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The  role  of  minerogenic  particles  in  light  scattering  in  several  lakes  and  a  river  (total  of  ten  sites)  in 
central  New  York,  which  represent  a  robust  range  of  scattering  conditions,  was  evaluated  based  on  an 
individual  particle  analysis  technique  of  scanning  electron  microscopy  interfaced  with  automated  x-ray 
microanalysis  and  image  analysis  (SAX),  in  situ  bulk  measurements  of  particle  scattering  and  backscat- 
tering  coefficients  ( bp  and  bhp ),  and  laboratory  analyses  of  common  indicators  of  scattering.  SAX  provided 
characterizations  of  the  elemental  x-ray  composition,  number  concentration,  particle  size  distribution 
(PSD),  shape,  and  projected  area  concentration  of  minerogenic  particles  (PAVm)  of  sizes  >0.4  jxm.  Mie 
theory  was  applied  to  calculate  the  minerogenic  components  of  bp  (6,„)  and  bbf)  (bbm)  with  SAX  data. 
Differences  in  PAV,„,  associated  primarily  with  clay  minerals  and  CaCO;J,  were  responsible  for  most  of  the 
measured  differences  in  both  bp  and  bhp  across  the  study  sites.  Contributions  of  the  specified  minerogenic 
particle  classes  to  bm  were  found  to  correspond  approximately  to  their  contributions  to  PAV,„.  The 
estimates  of  bln  represented  substantial  fractions  of  bp,  whereas  those  of  bb  m  were  the  dominant  compo¬ 
nent  of  bhp.  The  representativeness  of  the  estimates  of  bm  and  bbm  was  supported  by  their  consistency  with 
the  bulk  measurements.  Greater  uncertainty  prevails  for  the  bhm  estimates  than  those  for  b„n  associated 
primarily  with  reported  deviations  in  particle  shapes  from  sphericity.  The  PSDs  were  well  represented 
by  the  “£T  component  of  the  two-component  model  or  a  three  parameter  generalized  gamma  distribution 
[Deep-Sea  Res.  Part  I  40,  1459  ( 1993)1 .  The  widely  applied  Junge  (hyperbolic)  function  performed  poorly 
in  representing  the  PSDs  and  the  size  dependency  of  light  scattering  in  these  systems,  by  overrepre¬ 
senting  the  concentrations  of  submicrometer  particles  especially.  Submicrometer  particles  were  not 
important  contributors  to  bm  or  bb  m.  ©  2007  Optical  Society  of  America 
OCIS  codes:  120.5820,  180.5810,  290.1350,  290.4020,  290.5850. 


1.  Introduction 

Light  scattering  by  particles  is  a  fundamental  regu¬ 
lator  of  features  of  radiative  transfer  within  aquatic 
ecosystems  and  the  emergent  flux  signal  available  for 
remote  sensing,  and  it  is  important  in  assessments  of 
physical  properties  of  particles  [1].  The  magnitude 
and  character  of  particle  scattering  are  regulated  by 
four  features  of  a  particle  population:  particle  num¬ 
ber  concentration  (A0,  particle  size  distribution  (PSD) 
12-4],  particle  composition  [5],  and  particle  shape  [6]. 
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The  total  particle  scattering  coefficient,  bp(\)  (m  *), 
is  defined  as 


bp(\)  =  2tt  Pp(0,  X)sin  OdO,  (1) 


where  X  is  the  wavelength  of  the  light,  [3/;(0,  X)  is  the 
particle  volume  scattering  function  (VSF)  describing 
the  angular  distribution  of  the  scattered  radiation, 
and  0  is  the  scattering  angle;  integrating  p;,(0,  X) 
[Eq.  (1)]  over  the  angular  range  from  tt/2  to  tt  pro¬ 
vides  the  particle  backscattering  coefficient,  bhp(\). 
The  particle  backscattering  ratio  bhp(\)  corresponds 
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to  bbp(\)/bp(\).  The  bulk  (3p(9,  X),  fep(X),  and  bhp(\)  are 
all  inherent  optical  properties  (IOPs)  of  the  water 
medium  and  therefore  are  independent  of  the  geo¬ 
metric  structure  of  the  light  field. 

The  dependence  of  bp(\)  and  bbp{\)  on  particle  size 
fractions  and  compositional  components  remains  un¬ 
certain  for  marine  waters  [1,7,8],  while  related  infor¬ 
mation  for  inland  waters  is  extremely  rare  [9].  The 
prevailing  opinion  for  marine  systems  is  that  parti¬ 
cles  larger  than  1  |xm  in  size  contribute  the  most  to 
6P(X)  [1],  while  submicrometer  particles  (e.g.,  0.2 
-1  p,m)  regulate  bhp(\)  [1,8].  Moreover,  minerogenic 
particles  are  generally  believed  to  play  relatively 
more  important  roles  in  regulating  scattering  in  most 
case  2  marine  waters  [1,5,7],  and  probably  in  inland 
waters  [9],  compared  to  in  case  1  waters,  though  de¬ 
finitive  inorganic  particle  information  has  been  lack¬ 
ing  [8]. 

IOPs  are  additive;  the  total  is  the  sum  of  the  con¬ 
tributing  components  [10,11],  as  illustrated  here  for 
total  scattering,  6(X): 

6(\)  =  6P(X)  +  bw(k),  (2) 

where  btv(k)  and  6P(X)  are  the  scattering  coefficients  of 
pure  water  (seawater  for  marine  systems)  and  parti¬ 
cles,  respectively.  The  magnitude  of  bw(\)  is  small 
relative  to  bp(\)  in  the  vast  majority  of  waters  (i.e., 
b  ~  bp).  A  robust  representation  would  partition  bp(\) 
into  organic  [60(X)]  and  minerogenic  [6m(X)]  particle 
components  according  to 

6P(\)  =  60(X)  +  6m(X).  (3) 

Stramski  et  al.  [7,8]  supported  a  reductionist  approach 
of  increased  partitioning  of  IOPs  into  contributing 
components  to  develop  improved  understanding  of  the 
substantial  variability  in  ocean  optical  properties  and 
its  origins.  There  is  adequate  evidence,  though  more 
circumstantial  in  nature,  that  a  similar  approach 
would  be  necessary  to  understand  dynamics  within 
individual  freshwater  systems  [12],  as  well  as  differ¬ 
ences  among  systems  [13]  that  are  similar  to  varia¬ 
tions  among  case  2  systems  [5].  A  reductionist  version 
of  Eq.  (3)  would  be 

Nu  Nm 

MM  =  2  MM  +  2  M(M>  (4) 

1=1  7=1 

where  bo  i(\)  is  the  scattering  coefficient  of  the  ith  or¬ 
ganic  particulate  component  and  bmj(\)  is  the  scatter¬ 
ing  coefficient  of  the  yth  minerogenic  particulate 
component.  Application  of  the  reductionist  approach 
by  Stramski  et  al.  [7]  to  model  scattering,  by  using 
multiple  particulate  components  under  a  variety  of 
hypothetical  conditions,  yielded  insights  concerning 
optical  variability  in  the  ocean.  Their  analysis  focused 
on  partitioning  the  organic  particle  population,  i.e., 
fc0(X).  Twenty-one  components  were  included  in  the 
analysis:  18  planktonic  (one  for  viruses,  one  for  het- 
erotrophic  bacteria,  and  16  phytoplankton  groups/ 
species),  and  one  each  for  detritus,  mineral  particles, 
and  bubbles.  Their  speculative  model  simulations  de¬ 


picted  a  potentially  important  role  of  minerogenic 
particles  in  regulating  bp(\)  and  bbp(\)  in  both  case  1 
and  case  2  waters. 

Advancement  of  the  reductionist  approach  in  par¬ 
titioning  IOPs  has  been  limited  by  the  lack  of  tech¬ 
nologies  to  quantify  light-scattering  attributes  of 
multiple  particle  classes  in  real  ecosystems.  While 
measurements  of  N  and  PSDs  have  been  made 
widely  (e.g.,  Coulter  counters  [14]),  the  lack  of 
chemical  characterization  by  these  technologies 
prevents  partitioning  according  to  composition.  The 
size  threshold  of  many  of  these  instruments  (e.g., 
>1  pm)  is  also  problematic,  particularly  for  the  bbp(\) 
issue  [1,8].  Recently,  using  an  individual  particle 
analysis  technique,  scanning  electron  microscopy  in¬ 
terfaced  with  automated  x-ray  microanalysis  and  im¬ 
age  analysis  (SAX),  Peng  and  Effler  [15]  described  the 
light-scattering  attributes  of  minerogenic  particles 
(>  0.2  pm)  in  a  turbid  reservoir  where  these  particles 
dominate  6P(X).  SAX  provides  characterizations  of  el¬ 
emental  x-ray  composition,  N,  PSD,  and  shape  for 
multiple  geochemical  classes  of  minerogenic  parti¬ 
cles.  “Reasonably  good”  closure  was  reported  [15]  be¬ 
tween  Mie  theory  estimates  of  fe,„(66 0)  based  on  SAX 
data  and  measurements  of  the  beam  attenuation  co¬ 
efficient  at  660  nm  [c(660)  ~  1.045  X  6p(660)]  [5]. 

Here  we  describe  and  contrast  the  light-scattering 
features  of  several  lakes  and  a  river  based  on  indi¬ 
vidual  particle  analyses,  in  situ  bulk  measurements 
of  bp(\)  and  bhp(\),  and  laboratory  analyses  of  indica¬ 
tors  of  scattering.  We  use  SAX  to  characterize  the 
light-scattering  attributes  of  minerogenic  particles  in 
these  systems  and  partition  this  fraction  into  note¬ 
worthy  contributing  generic  type  components.  Imple¬ 
mentations  of  the  forward  method  in  modeling  light 
scattering  and  of  the  reductionist  approach  in  parti¬ 
tioning  IOPs  are  advanced  through  Mie  theory  cal¬ 
culations  of  bm{\)  and  &6m(X).  The  findings  are  used  to 
evaluate  and  contrast  the  role  of  minerogenic  compo¬ 
nents  in  regulating  features  of  scattering  for  the 
study  systems  and  to  estimate  6„(X). 

2.  Methods 

A.  Study  System 

The  study  systems  include  five  lakes  and  a  21  km 
reach  of  a  river  located  in  central  New  York  (Fig.  1). 
These  alkaline  hard-water  systems  represent  a  sub¬ 
stantial  range  of  trophic  state  and  sediment  content. 
All  the  study  lakes  stratify  in  summer,  though  wide 
differences  in  flushing  rates  prevail  among  these  sys¬ 
tems  (Table  1).  Owasco,  Skaneateles,  and  Otisco 
Lakes  are  the  easternmost  of  New  York’s  11  Finger 
Lakes.  Onondaga  Lake,  located  in  metropolitan  Syr¬ 
acuse,  is  polluted  by  municipal  effluent  and  residual 
industrial  waste  [13].  Pelagic  lake  sites  were  as¬ 
sessed  (Fig.  1),  except  for  a  shallow  (~3  m)  unstrati¬ 
fied  location  in  the  southern  end  of  Otisco  Lake  that 
adjoins  the  mouth  of  the  largest  tributary  and  is  iso¬ 
lated  from  the  main  portion  of  the  lake  by  a  con¬ 
structed  causeway. 
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Fig.  1.  Study  area  and  sites  in  central  New  York.  See  Table  1  for 
site  abbreviations. 


All  the  study  lakes  drain  into  the  Seneca  River 
(basin  area  of  9000  km2,  average  flow  of  96  m3  s'1), 
which,  after  combining  with  the  Oneida  River  to  form 
the  Oswego  River  (Fig.  1),  enters  Lake  Ontario.  Large 
sediment  and  nutrient  loads  are  carried  by  the  river 
upstream  of  Cross  Lake  [16,17].  Cross  Lake  is  an 
intervening  lake  for  this  river  (>  95%  of  inflow  to  the 
lake)  [16].  Short-circuited  flow  of  the  river  input  to 
this  lake’s  outflow  (i.e.,  downstream  portions  of  the 
river),  which  promotes  differences  in  various  charac¬ 
teristics  along  the  lake’s  primary  axis  [16],  is  com¬ 
mon.  Two  sites  were  assessed  in  Cross  Lake  (Fig.  1), 
one  in  the  northern  portion,  the  second  further  south 


and  closer  to  the  river  inflow.  The  three  river  sites 
(adjoining  navigation  buoys  B424,  B409,  and  B317, 
Fig.  1)  were  selected  to  depict  the  effects  of  Cross 
Lake  and  Dreissena  polymorpha  (zebra  mussels)  me¬ 
tabolism,  respectively.  Dense  populations  of  this  bi¬ 
valve  invader  prevail  between  Cross  Lake’s  outflow 
and  the  downstream  study  boundary  (between  buoys 
B409  and  B317),  resulting  in  substantial  changes  in 
common  measures  of  water  quality  over  this  reach, 
including  increased  clarity  [17]. 

B.  Sampling,  and  Laboratory  and  Field  Measurements 

Near  surface  samples  were  collected  with  a  submers¬ 
ible  pump  between  1000  and  1400  h  on  four  days  in 
July  2005.  Common  indicators  of  scattering  and  wa¬ 
ter  quality  were  measured  in  the  laboratory,  includ¬ 
ing  c(660)  (10  cm  path  length  transmissometer; 
C-Star/WET  Labs),  turbidity  ( T,n  nepholometric  tur¬ 
bidity  units  (NTUs);  Hach  Model  2100AN;  Clesceri 
et  al.  [18],  total  and  fixed  suspended  solids  (TSS  and 
FSS;  1.5-p.m  pore  sized  Whatman  glass  fiber  filter  for 
TSS,  and  FSS  after  550  °C,  i.e.,  inorganic),  and  vol¬ 
atile  suspended  solids  (VSS  (=TSS  -  FSS),  approxi¬ 
mately  organic  [18]),  and  chlorophyll  a  (Chi;  Parsons 
et  al.  [19]).  TSS  is  equivalent  to  suspended  particu¬ 
late  matter  (SPM),  the  term  commonly  adopted  in 
marine  studies.  To  account  for  the  underestimation  of 
FSS  due  to  the  structural  water  loss  in  clay  minerals 
during  the  ignition  process,  a  fraction  (structural  wa¬ 
ter  mass  in  argillaceous  minerals)  of  0.09  [20]  was 
used  to  estimate  the  water  content  in  the  clay-FSS 
(based  on  contributions  of  projected  area  concentra¬ 
tions  of  clay  particles  to  that  of  total  mineral  parti¬ 
cles;  see  SAX  protocols).  This  water  loss  was  then 
subtracted  from  and  added  to  VSS  and  FSS,  respec¬ 
tively,  to  correct  the  gravimetrically  measured  parti¬ 
cle  concentrations  (VSScorr  and  FSScorr). 

Tn  and  c(660)  are  surrogate  measures  of  bp  [5,10]. 
Measurements  of  c(660)  and  Tn ,  and  filtering  for  TSS 


Table  1.  Characteristics  of  Study  Systems  and  Location  of  Sites 


System" 

Location6 

Mean 

Depth 

Surface 

Area 

Trophic 

Sampling 

Date 

(Abbreviation) 

Latitude 

Longitude 

(m) 

(km2) 

State' 

(July) 

Owasco  L.  (Ow) 

42°  52'03.35"  N 

76°  31'24.50"  W 

29.3 

26.7 

m 

28 

Skaneateles  L.  (Sk) 

42°  52'05.35"  N 

76°  31 '24.50"  W 

43.5 

35.9 

0 

8 

Otisco  L. 

10.2 

7.6 

m 

North  (OtN) 

42°  52 '2 1.24"  N 

76°  17'42.10"  W 

— 

— 

— 

28 

South  (OtS) 

42°  50'40.77"  N 

76°  15'22.33"  W 

— 

— 

— 

28 

Onondaga  L.  (On) 

43°  06 '54"  N 

76°  14 '34"  W 

10.9 

12.0 

e 

27 

Cross  L. 

5.5 

9.0 

e 

North  (CrN) 

43°  08'23.47"  N 

76°  29 '05. 94"  W 

— 

— 

— 

29 

South  (CrS) 

Seneca  R. 

43°  06  57.01"  N 

76°  28'28.67"  W 

— 

— 

— 

29 

B424 

43°  06'05.02"  N 

76°  29'57.21"  W 

— 

— 

— 

29 

B409 

43°  06' 14. 18"  N 

76°  26'41.19"  W 

— 

— 

— 

29 

B317 

43°  08'49.17"  N 

76°  18'50.16"  W 

— 

29 

r*See  Fig.  1. 

'’World  Geodetic  System  1984. 

‘o,  oligotrophic;  m,  mesotrophic;  e,  eutrophic. 
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and  Chi,  were  conducted  on  the  day  of  sample  collec¬ 
tion.  Measurements  of  c(660)  were  also  made  in  situ 
(same  model  as  laboratory  unit)  at  approximately  the 
time  of  sampling  with  profiling  instrumentation 
(SeaBird  SeaLogger  Profiler).  Additionally,  spectral 
absorption,  attenuation,  and  particulate  backscatter- 
ing  coefficients  [a(X),  c(X),  and  bbp(\)]  were  measured 
using  a  combined  instrument  package  (same  profiling 
cage,  integrated  outputs)  of  ac-s  (spectral  absorption 
and  attenuation  meter,  path  length  25  cm)  and  BB9 
(WET  Labs,  Inc.,  Philomath,  Oregon,  USA),  also  at 
approximately  the  time  of  sampling.  The  ac-s  mea¬ 
sures  a(X)  and  c(X)  (relative  to  pure  water),  and  has 
a  spectral  resolution  of  4  nm  over  the  range  of 
400-730  nm.  The  instrument  was  calibrated  accord¬ 
ing  to  guidelines  of  the  manufacturer  [21]  with  pure 
water  and  checked  between  deployments  by  an  air 
calibration  protocol.  Corrections  for  differences  in 
pure-water  absorption  and  attenuation  due  to  tem¬ 
perature  were  made  according  to  Pegau  et  al.  [22]. 
The  ac-s  absorption  spectra  were  corrected  for  scat¬ 
tering  error  through  a  WET  Labs  protocol  that  uses  a 
reference  X  to  determine  the  value  of  the  scattering 
coefficient  to  be  subtracted  from  a(X).  The  spectral 
particulate  scattering  coefficient,  6P(X),  was  obtained 
as  the  difference,  c(X)  -  a(X).  The  relative  errors  for 
attenuation  measurements  as  a  result  of  multiple 
scattering  (corresponding  to  optical  thickness  of  <0.8 
for  the  C-Star  and  <2  for  the  ac-s)  were  <5%  [23], 
with  bhp{\)  obtained  at  nine  wavelengths  (412,  440, 
488,  510,  532,  595,  650,  676,  and  715  nm)  by  the  BB9 
backscattering  meter.  This  instrument  measures  the 
VSF  at  the  angle  of  117°  [p(117°,  X)],  an  angle  iden¬ 
tified  as  a  minimum  convergence  point  for  variations 
in  p(0,  X)  caused  by  suspended  particles  and  water 
itself.  Corrections  for  attenuation  along  the  path  from 
the  light  source  to  the  sample  volume  to  the  detector 
were  based  on  the  paired  ac-s  measurements  of  a(X), 
using  a  protocol  described  by  the  manufacturer  [24]. 
Particulate  VSF  was  measured  at  117°,  (3^(117°,  X) 
=  [3(117°,  X)  -  (3U,(117°,  X),  the  latter  being  the  vol¬ 


ume  scattering  of  water.  The  value  of  bbp(\)  was  es¬ 
timated  from  the  (3^(117°,  X)  according  to,  bbp(\)  = 
2'ttxPp(117°,  X),  where  the  value  of  x  is  1.1  [25]. 

C.  Scanning  Electron  Microscopy  Interfaced  with 
Automated  X-Ray  Microanalysis  and  Image  Analysis 
Protocols 

Sample  preparation  and  analytical  methods  for  SAX 
have  been  described  previously  [26,27];  salient  fea¬ 
tures  are  reviewed  here.  Particles  were  deposited 
onto  polycarbonate  membranes  (25  mm  diameter, 
0.4  pm  pore  size)  by  pressure  filtration  and  carbon 
coated.  SAX  was  conducted  with  an  Aspex  Microanal¬ 
ysis  System  with  Perception  Analysis  Software  (with 
a  plug-in  Automated  Feature  Analysis  Module).  Ap¬ 
proximately  2000  particles  were  analyzed  for  each  of 
the  ten  samples  for  both  elemental  and  morphometric 
characteristics.  The  time  required  for  each  analysis 
was  ~3  h.  Particles  of  sizes  >~0.2  ^m  are  character¬ 
ized  (based  on  the  electron  beam  search  step  size),  but 
only  those  greater  than  the  nominal  filter  pore  size 
are  reported  here. 

Composition  of  individual  particles  is  assessed  by 
SAX  based  on  the  acquisition  of  x-ray  counts  for  16 
elements  (Na  and  higher  atomic  number,  including 
Al,  Si,  Ca,  Fe,  Mn).  Elemental  composition  of  a  par¬ 
ticle  is  represented  by  its  normalized  elemental  x-ray 
counts,  i.e.,  relative  intensities  [28] .  Particles  were 
partitioned  into  eight  generic  biogeochemical  types  or 
classes  (Table  2)  according  to  a  scheme  developed 
earlier  for  the  Finger  Lakes  of  New  York  [27];  the 
definitions  of  particle  types  are  consistent  with  the 
geochemical  and  mineralogical  setting  of  the  region. 
This  classification  scheme  differs  from  that  used  for 
the  low  alkalinity  (and  low  [Ca2*])  waters  of  the 
Catskill  region  [26]  of  New  York  by  the  inclusion  of 
types  for  particles  containing  Ca,  consistent  with  the 
summertime  precipitation  of  CaC03  (calcite)  that  is 
known  to  occur  in  several  of  the  lakes  of  this  study 
[13].  The  “Clay”  class  corresponds  to  clay  minerals, 


Table  2.  Specification  of  Generic  Particle  Types,  According  to  X-Ray  Characteristics,  and  the  Complex  Refractive  Indices  (Relative  to  Water)  Used  in 

Mie  Calculations 


Type 

Description 

X-Ray  Characteristics 

Sources/Origins 

m 

Organics 

Biological 

Low  x-ray  net  counts  (<750)“ 

Autochthonous/terrigenous 

Clay 

Aluminosilicates 

Al,  5  to  55%;  Si,  20  to  85%; 

Al  plus  Si  >  50%h 

Terrigenous 

1.173-0.001/ 

Ca-rich 

CaC03  precipitates 

Ca  >  85% 

Autochthonous 

1.199-0.0001/ 

Ca-agg 

CaC03  coating  on 
other  particles 

20%  <  Ca  <  85% 

Autochthonous/terrigenous 

1.186-0.0005/ 

Quartz 

Mineral  Si02 

Si  >  90%,  high  x-ray  density* 

Terrigenous 

1.155-0.0001/ 

Diatom 

Biogenic  Si0.2 

Si  >  90%,  low  x-ray  density* 

Autochthonous 

Si-rich 

Si-containing 
minerals,  silicates 

60%  <  Si  <  90% 

Terrigenous 

1.173-0.001/ 

Other*' 

Miscellaneous 

particles 

Not  specified 

Various 

1.173-0.001/ 

"Live  x-ray  acquisition  time  was  3  s. 

'’All  percentages  are  elemental  x-ray  relative  intensities. 

‘  The  x-ray  density  refers  to  the  ratio  of  a  particle’s  total  x-ray  counts  to  its  size  [26]. 
‘'This  incorporates  all  inorganic  particles  not  captured  in  the  specified  classes. 
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“Quartz”  to  Si02  (mineral  phase),  and  “Si-rich”  to 
other  silicate  minerals.  These  classes  have  terrig¬ 
enous  origins.  The  “Diatom”  and  “Quartz”  particles 
were  differentiated  according  to  the  protocol  of  Peng 
et  al.  [26],  which  is  based  on  the  lower  x-ray  count 
rates  and  the  larger  size  of  diatoms  compared  to  those 
of  quartz.  The  “Ca-rich”  class  corresponds  to  calcite 
particles  or  CaC03-coated  particles,  and  reflects  au¬ 
tochthonous  formation  in  the  water  columns  of  the 
study  systems  [29].  The  “Ca-agg”  group  represents 
aggregate  particles  that  have  been  attributed  to  the 
coating  of  terrigenous  particles  by  calcite  precipita¬ 
tion  [28].  The  “Other”  class  contains  inorganic  parti¬ 
cles  not  represented  by  the  specified  groups  (Table  2). 

The  “Organics”  class  systematically  underrepre¬ 
sents  the  contributions  of  organic  particles  due  to  their 
common  composition  of  low  atomic  number  elements 
(lower  than  that  of  Na),  as  a  result  of  methodological 
limitations  inherent  to  SAX  [26,28].  Nonetheless,  the 
enumeration  of  this  particle  population  by  SAX  and 
subsequent  calculations  provide  semi-quantitative  es¬ 
timates  of  its  contribution  to  the  overall  PAV. 

Morphometric  characterization  of  particles  by  SAX 
is  based  on  a  “rotating  chord”  algorithm,  which  pro¬ 
vides  a  detailed  representation  of  morphometries  of 
individual  particles  [15],  including  projected  area 
(PA),  length,  and  width.  Particle  size  (< d )  is  calculated 
from  PA  as  the  circular  area  equivalent  diameter. 
The  total  projected  area  concentration  (total  sample 
PA  per  unit  volume  of  water,  PAV)  is  the  sum  of  the 
individual  PA  values  (result  of  extrapolation  based 
on  the  fraction  of  analyzed  filter  area)  divided  by 
the  filtered  volume  of  water  sample.  PAV  can  be 
partitioned  into  contributing  components,  e.g., 
minerogenic  PAV  (PAVm),  clay-PAV,  or  other  type 
components  as  listed  in  Table  2. 

The  “nonsphericity”  of  a  particle  (e.g.,  shape  factor) 
is  represented  by  the  aspect  ratio  (ASP),  defined  as 
the  ratio  of  a  particle’s  length  to  its  width  [30].  A 
sphere  has  an  ASP  value  of  1;  increases  in  ASP  cor¬ 
respond  to  greater  deviations  from  spheres.  Morpho¬ 
metric  features  were  used  to  resolve  the  contribution 
of  Phacotus  (a  phytoplankton  with  a  CaC03  lorica 
[27])  to  the  Ca-rich  class,  based  on  its  nearly  spher¬ 
ical  shape  (i.e.,  ASP  approaches  1)  and  large  size 
(d  >  8  pm).  This  differentiation  effectively  results  in 
nine  delineated  particle  classes. 

D.  Particle  Size  Distributions,  Calculations  of  bm  and 
bb  m,  and  Estimates  of  b0 

PSDs  are  presented  as  F(d),  such  that  F{d)i\d  gives 
the  PAV  of  water  in  the  size  range  d  ±  l/2(AcD  (d  is 
the  midpoint  and  Ad  is  the  width  of  a  size  bin).  The 
measured  minerogenic  PSDs  were  fitted  with  two 
different  functional  representations,  the  Junge  (or 
hyperbolic)  and  the  B  component  of  the  two- 
component  (2C-R)  model.  The  fitting  of  the  Junge 
model  [F(d)  =  Cd  J ,  where  C  is  a  constant  depending 
on  N ,  and  j  is  the  slope  of  the  distribution  when 
linearity  is  assumed  in  the  log-log  plot  of  F(d)],  was 
limited  to  particles  with  d  >  1.25  pm,  a  typical  lower 


threshold  for  particle  counters  [2,31].  Risovic’s  two- 
component  model  (TCM)  [4,32]  is  defined  as 

F(d)  =  CAd»A  exp(-52 dyA)  +  C exp(- 17 tf™),  (5) 

where  CA  and  CB  are  proportional  to  the  number 
concentrations  of  small  and  large  particles,  respec¬ 
tively,  and  p,  and  y,  (i  =  A,  B)  are  parameters  of  the 
distribution.  In  our  earlier  application  of  this  model, 
we  have  found  that  the  B  component  alone  ade¬ 
quately  represented  the  particle  sizes  (>0.5  pm)  re¬ 
sponsible  for  6m(660)  in  Schoharie  Reservoir,  New 
York,  USA  [15];  this  reduced  the  number  of  function 
coefficients  by  half,  and  avoided  negative  values  of  CA 
when  Eq.  (5)  was  fitted  to  SAX-measured  PSD  data 
from  the  reservoir  [15]. 

The  bulk-scattering  coefficients  for  minerogenic 
particles  were  calculated,  based  on  the  individual 
particle  information  from  SAX,  according  to 

Nm 

b,n(K)  =  2  QbmAmi>  x-  (6) 

J  =  1 

where  Nrn  is  the  number  of  minerogenic  particles  per 
unit  volume  of  water  (number  m  '*),  Q/)m,,  is  the  di¬ 
mensionless  efficiency  factor  for  scattering  of  the  ith 
minerogenic  particle  with  the  adjusted  projected  area 
of  PA*  (  (m2;  see  subsequent  text);  likewise,  bhm(\) 
was  calculated  with  the  backscattering  efficiency  fac¬ 
tor,  Qbbm,r  The  values  of  Q  factors  depend  on  the 
particle’s  relative  (to  water)  refractive  index  (m  =  n 
-  in',  where  n  and  n'  are  the  real  and  imaginary 
parts  of  the  complex  index  of  refraction,  respectively), 
its  size  (d,),  and  \.  The  results  are  reported  for  \ 
=  650  nm  [i.e.,  6„J(650),  6fem(650)],  a  wavelength  for 
which  bh{\)  measurements  were  made,  and  close  to 
the  reference  wavelength  (660  nm)  adopted  by  Peng 
and  Effler  [15]  in  similar  calculations.  Nonminero- 
genic  particles  (i.e.,  those  of  the  Organics  and  Diatom 
classes,  and  the  Phacotus  fraction  of  the  Ca-rich 
class)  were  excluded  in  these  scattering  calculations. 

We  calculated  Qhmtl  and  Qbbmj  for  individual  min¬ 
eral  particles  using  the  Mie  theory  algorithm  (for 
homogeneous  spheres)  provided  by  Bohren  and  Huff¬ 
man  [33].  Values  of  n  (Table  2)  were  specified  or 
guided  by  the  listings  presented  by  Kerr  and  Rogers 
[34]  and  Wozniak  and  Stramski  [35].  The  value  for 
the  Ca-rich  class  corresponded  to  that  of  calcite.  The 
Clay,  Si-rich,  and  Other  classes  were  assigned  the 
value  of  kaolinite,  whereas  the  Ca-agg  class  corre¬ 
sponded  to  the  average  of  kaolinite  and  calcite.  The  n' 
values  assigned  to  the  particle  types  (Table  2)  were 
consistent  with  the  bounds  adopted  in  analyses  by 
Wozniak  and  Stramski  [35],  and  reflected  the  as¬ 
sumption  of  minimal  absorption  by  Ca-rich  and 
Quartz  classes  compared  to  the  others.  The  calcu¬ 
lated  contributions  of  various  size  classes  of  minero¬ 
genic  particles  to  the  overall  estimated  scattering  are 
presented  in  a  cumulative  format  (e.g.,  [5]).  Patterns 
obtained  with  functional  fits  are  compared  to  those 
based  on  observations.  Sensitivity  analyses  were  con¬ 
ducted  (for  observed  d  range  of  each  sample)  to  depict 
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the  dependence  of  scattering  on  n  values  for  three 
PSD  cases:  the  measured,  and  two  described  by  the 
Junge  function;  values  of  6m(650)  and  bb m(650)  were 
calculated  for  these  cases  of  PSD  with  uniform  values 
of  n  from  1.05  to  1.25  in  increments  of  0.05.  The  two 
Junge  function  scenarios  corresponded  first  to  the 
fitted  j  (for  observations  d  >  1.25  |jim)  and  second 
specified  j  =  4. 

Adjustments  in  PA  values  determined  for  non- 
spherical  particles  by  SAX  are  needed  to  represent 
their  effective  projected  area  (PA*)  within  the  water 
column,  because  of  their  nonrandom  orientation 
(“lying  flat”)  on  the  filter  [14,15].  A  PA*/PA  ratio 
value  of  0.87  was  reported  for  quartz  dust  particles  by 
Proctor  and  Harris  [36],  and  was  adopted  in  an 
analysis  of  light  scattering  by  Jonasz  [14].  The  light¬ 
scattering  closure  analysis  based  on  SAX  observa¬ 
tions  for  Schoharie  Reservoir,  where  clay  minerals 
dominated,  supported  a  PA*/PA  value  of  —0.7  for 
those  particles  [15].  The  effect  of  overestimation  of 
PA  on  6m(650)  and  6fom(650)  was  bracketed  here  by 
two  sets  of  adjustments:  first  for  particles  with  ASP 
values  >1.5,  then  that  of  >2.5,  adopting  a  PA*/PA 
ratio  of  0.8.  The  two  ASP  thresholds  were  selected  to 
reflect  reasonable  upper  and  lower  bound  adjust¬ 
ments  for  the  effect.  Corresponding  adjustments  in 
d(d *)  were  made  in  these  Mie  theory  calculations. 

The  organic  scattering  component,  fc„(6 50),  was  cal¬ 
culated  as  the  difference  between  the  in  situ  mea¬ 
surements  of  bp(6 50)  and  the  SAX-based  estimate  of 
6m(650).  Additionally,  independent  estimates  of  b() 
were  based  instead  of  Chi  concentration  (6o  Chl)  using 
a  relationship  developed  for  case  1  waters  [37],  mod¬ 
ified  to  represent  an  inverse  dependency  on  \  [38], 

60iChi(X)  =  0.347[Chl]0766(660/\)0.97.  (7) 

The  average  scattering  and  backscattering  efficien¬ 
cies  of  the  minerogenic  particle  populations,  (Q6>m) 
and  ( Qbb,m)>  were  calculated  as  the  ratios  of  the  es¬ 
timates  of  minerogenic  scattering  (6m  and  bhm)  to 
PAVm  (as  the  sum  of  measured  PA,,  without  adjust¬ 
ment).  In  an  analogous  manner,  average  efficiencies 
were  also  calculated  for  the  various  contributing  par¬ 
ticle  classes. 

3.  Results 

A.  Indicators  of  Light  Scattering  and  Contributing 
Constituents 

Wide  differences  in  light-scattering  levels  were  de¬ 
picted  by  the  common  metrics  of  water  quality  for  the 
various  study  sites  (Fig.  2).  Values  of  c(660)  and  Tn 
ranged  from  0.48  m  1  and  0.4  NTU  for  Skaneateles 
Lake  to  14.7  m  1  and  22  NTU  for  the  south  site  on 
Otisco  Lake  [Fig.  2(a)].  These  two  surrogate  metrics 
of  scattering  were  strongly  correlated  (r  =  0.97, 
p  <  0.01).  Concentrations  of  TSS  ranged  from  —0.3 
(Skaneateles  Lake)  to  nearly  20  g  m  3  [south  site  on 
Otisco  Lake;  Fig.  2(b)].  Wide  differences  in  the  frac¬ 
tion  FSScorr/TSS  prevailed,  ranging  from  0.34  for  On¬ 
ondaga  Lake  to  0.86  for  B424  on  the  Seneca  River; 


Ow  Sk  OtN  OtS  On  B424  CrN  CrS  B409B317 


Fig.  2.  Bulk  water  quality  metrics  for  study  sites  related  to  light 
scattering:  (a)  c(660)  and  turbidity,  (b)  suspended  solids,  and  (c) 
chlorophyll  concentration,  Chi.  See  site  abbreviations  in  Table  1. 

inorganic  components  also  made  large  contributions 
to  TSS  at  the  B409  river  site  and  the  south  site  on 
Otisco  Lake.  Chi  values  ranged  from  0.7  (Skaneateles 
Lake)  to  nearly  50  mg  m  5  [north  site  on  Cross  Lake 
and  B409  on  the  river;  Fig.  2(c)].  Concentrations  ex¬ 
ceeded  10  mg  m  3  (i.e.,  bloom  conditions)  for  five  of 
the  ten  sites. 

Empirical  analyses  established  that  TSS  was  a 
good  predictor  of  differences  in  light  scattering 
among  the  study  sites;  the  linear  least-squares  re¬ 
gression  relationship, 

c(660)  =  (0.713  x  TSS)  +  0.613,  (8) 

explained  96%  (p  <  0.001)  of  the  differences  in  c(660). 
A  similar  relationship  explained  98%  (p  <  0.001)  of 
the  differences  in  Tn.  Relationships  to  explain  the 
differences  in  these  surrogate  metrics  of  light  scat¬ 
tering  based  on  Chi  were  weak  and  insignificant. 
However,  a  relationship  that  partitioned  the  two 
components  of  TSS, 

c(660)  =  (0.558  x  FSScorr)  +  (1.227  x  VSScorr)  +  0.105, 

(9) 

performed  equally  well  (r2  =  0.98)  as  the  single  com¬ 
ponent  representation.  The  coefficients  for  both 
FSScoit  and  VSScorr  were  significant  (p  <  0.001),  but 
the  intercept  was  not.  Variations  in  FSScorr  explained 
88%  of  the  overall  variability  in  c(660);  the  addition  of 
VSScorr  explained  another  10%.  This  suggests  that 
differences  in  the  observed  light  scattering  were  pri- 
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marily  due  to  minerogenic  particle  populations  and 
secondarily  to  organic  particulate  components. 

B.  Composition  and  Morphometries  of  Particle 
Populations 

The  composition  of  particle  populations  is  presented 
in  the  context  of  contributions  of  particle  type  com¬ 
ponents  to  PAV  (Table  3)  because  of  the  central  role 
this  attribute  plays  in  light  scattering.  The  adopted 
combined  geochemical  and  physical  classification 
scheme  performed  well  in  representing  the  minero¬ 
genic  particles  of  these  systems,  as  less  than  6%  of  the 
PAV  was  associated  with  the  Other  category.  PAVm 
measurements  were  strongly  correlated  (r  =  0.97) 
with  FSScorr  observations.  Clay  minerals  [Fig.  3(a)], 
CaC03  [Figs.  3(b)  and  3(d)],  and  CaC03  aggregates 
[Fig.  3(a)],  common  to  all  the  study  systems,  were  the 
most  important  minerogenic  components.  The  Clay 
class  was  the  single  largest  component  of  PAVm  at 
seven  of  the  sites,  whereas  CaC03  (non -Phacotus 
component  of  Ca-rich,  together  with  Ca-agg)  was  the 
largest  component  for  the  others.  Both  the  Quartz 
[Fig.  3(a)]  and  Si-rich  classes  made  smaller  contribu¬ 
tions.  The  Clay  class  made  the  greatest  contribution 
at  B424,  upstream  of  Cross  Lake,  representing  68.5% 
and  71.7%  of  the  total  and  minerogenic  PAV,  respec¬ 
tively  (Table  3).  This  class  made  larger  contributions 
to  PAV  at  the  south  site  on  Cross  Lake,  which  was 
closer  to  the  river  inflow,  than  at  the  north  site.  The 
clay  mineral  PAV  and  its  contribution  were  much 
larger  at  the  south  shallow  site  on  Otisco  Lake  than 
for  the  north  pelagic  site.  Calcium-containing  particle 
types  were  the  dominant  contributors  to  PAV  at  the 


northern  site  of  Otisco  Lake  (Ca-rich  plus  Ca-agg, 
—58%  of  PAV;  65%  of  PAVJ  and  Owasco  Lake.  Di¬ 
atoms  dominated  PAV  in  Skaneateles  Lake  and  were 
a  noteworthy  contributor  in  Cross  Lake  and  in  down¬ 
stream  portions  of  the  river.  Contributions  by  Pha¬ 
cotus  were  noteworthy  (e.g.,  >6%)  in  Cross  Lake  and 
the  south  site  on  Otisco  Lake.  The  conspicuously 
higher  PAV  values  for  the  Organics  class  for  the 
north  Cross  Lake  and  Onondaga  Lake  sites  sug¬ 
gested  greater  contributions  to  scattering  by  organic 
particles  at  these  locations. 

The  ASP  values  for  the  Diatom  class  and  Phacotus 
fraction  of  the  Ca-rich  class  present  a  valuable  qual¬ 
itative  context  for  considering  the  shape  features  of 
the  minerogenic  particles  (Table  4).  Phacotus  [Fig. 
3(d)]  had  a  nearly  spherical  shape,  while  the  diatoms 
deviated  the  most  from  spherical,  associated  primar¬ 
ily  with  pennate  forms  such  as  Asterionella  sp.  [Fig. 
3(c)].  The  differences  in  average  ASP  values  for  the 
diatom  populations  suggest  differences  in  this  com¬ 
munity’s  composition  among  the  study  sites.  Average 
ASP  values  for  minerogenic  particle  classes  for  all 
sites  were  greater  than  one,  depicting  deviations 
from  spherical  shapes  (e.g.,  Fig.  3).  However,  these 
did  not  approach  the  much  higher  values  of  diatoms. 
Substantial  variability  in  shape  within  each  minero¬ 
genic  class  was  observed  at  each  site,  as  represented 
by  standard  deviations  about  the  mean  ASP  values 
(Table  4),  but  this  did  not  depend  on  size.  The  ASP 
values  for  the  Clay  and  Ca-agg  classes  were  some¬ 
what  higher  than  for  the  other  minerogenic  classes  at 
most  of  the  sites. 


Table  3.  Summary  of  PAV  Measurements  for  Study  Sites  According  to  Geochemical  Classes  with  Contributions  to  Total  PAV  and  PAV„ 


Type  Contribution 

(%) 

Low 

X-Ray 

Counts 

Organics 

Biogenic/Inorganic 

Coating 

Minerogenic 

PAV 
(cm2  1 

PAV,,, 

*)  (cm2 1'1) 

Ca-Rich 

System" 

Organics 

Diatom 

Phacotus 

Non -Phacotus 

Ca-Agg 

Clay 

Quartz 

Si-Rich 

Other 

Seneca  R. 

B424 

32.3 

30.8  (95.5/’ 

3.2 

1.4 

0.0 

3.6  (3.8)1 

13.6  04.2) 

68.5(71.7) 

4.4  (4.6) 

3.2  (3.4) 

2.2  (2.3) 

CrN 

10.8 

5.1(47.4) 

25.4 

15.5 

11.6 

15.3(32.3) 

13.1(27.7) 

9.9  (20.9) 

1.1  (2.3) 

4.5  (9.5) 

3.5  (7.4) 

CrS 

12.3 

8.7(70.3) 

6.9 

16.6 

6.4 

8.5(12.1) 

13.6  09.3) 

35.7  (50.8) 

2.0  (2.8) 

7.6  00.8) 

2.9(40) 

B409 

23.2 

17.7(76.6) 

8.1 

14.8 

1.2 

3.7  (4.8) 

11.2  04.6) 

47.2(61.6) 

4.8  (6.3) 

5.0  (6.5) 

4.7(60) 

B317 

3.1 

2.4  (77.9) 

4.1 

18.0 

0.0 

4.8  (6.2) 

6.8  (8.7) 

50.5  (64.8) 

9.6  02.3) 

2.2  (2.8) 

4.0(50) 

Lakes 

Ow 

2.4 

1.9(78.3) 

12.2 

6.6 

3.0 

27.3(34.8) 

14.108.0)  21.3(27.1) 

3.6  (4.6) 

6.4  (8.2) 

5.6(70) 

Sk 

2.0 

0.43  (21.4) 

0.6 

78.0 

0.0 

0.75  (3.5) 

2.0  (9.1) 

12.4  (58.0) 

1.8  (8.6) 

1.6  (9.5) 

2.8  03.3) 

OtN 

7.8 

7.0  (89.1) 

10.2 

0.1 

0.6 

28.9  (32.4) 

29.4  (33.0) 

25.7  (28.8) 

1.2  0.3) 

1.1  (1.2) 

2.8  (7.4) 

OtS 

50.9 

38.3  (75.2) 

5.4 

10.5 

8.8 

2.7  (3.6) 

5.5  (6.2) 

46.7(62.1) 

5.9  (7.8) 

8.8  01.7) 

5.61(7.4) 

On 

12.3 

7.1(57.7) 

41.8 

0.5 

0.0 

6.5(11.3) 

13.1  (22.7) 

28.8  (49.9) 

2.5  (4.3) 

1.8(30) 

5.0  (8.7) 

"See  Table  1  for  system  abbreviations. 

'’Percentages  of  PAVm  in  PAV  are  listed  in  parentheses. 

'  Numbers  in  parentheses  are  the  minerogenic  type  percentages  in  PAV,,,. 
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Fig.  3.  Micrographs  of  particles  from  study  sites:  (a)  Seneca 
River,  quartz  (cf  —  2  p.m)  in  upper  left,  clay  mineral  particles  ( d 
~  6-8  pm)  lower  right  and  upper  third/left  of  center,  and  Ca-agg 
particle  (d  ~  2  pm)  lower  half/left  of  center,  (b)  Onondaga  Lake, 
large  (>15  pm)  centric  diatom  (on  edge)  on  the  left,  calcite  particle 
(d  ~  8  pm)  on  lower  right,  (c)  Onondaga  Lake,  algae,  two  centric 
diatoms,  a  pennate  diatom  ( Asterionella  sp.)  in  center,  and 
Scenedesmus  sp.  (Organics)  center/right,  and  (d)  Owasco  Lake, 
Phacotus  sp.  (d  ~  10  pm)  lower  left,  dendritic  calcite  particle  upper 
center,  and  large  filamentous  diatom  (centric)  on  right. 


C.  Particle  Size  Distributions  of  Mineral  Particles 

The  minerogenic  PSDs  are  presented  in  two  groupings: 
lake  sites  [excluding  Cross  Lake,  Fig.  4(a)],  and  the 
Cross  Lake  and  Seneca  River  system  [Fig.  4(b)].  The 
highest  particle  concentrations  occurred  at  sizes  be¬ 
tween  0.4  and  0.5  |xm,  close  to  the  lowest  bound  of 
quantification,  with  the  conspicuous  exception  of  the 
north  site  on  Otisco  Lake,  where  the  maximum  was 


at  a  size  of  —  1.2  jxm  [Fig.  4(a)].  Progressive  decreases 
in  concentrations  were  observed  for  increased  parti¬ 
cle  sizes.  Much  higher  minerogenic  N  values  pre¬ 
vailed  throughout  the  0.4-10  |xm  size  range  at  the 
south  Otisco  Lake  site,  compared  to  the  pelagic  site  in 
the  northern  portion  of  the  lake.  Similar  PSD  pat¬ 
terns  were  observed  for  Onondaga,  Owasco,  and  Ska- 
neateles  Lakes,  though  substantially  higher  N  levels 
prevailed  in  Onondaga  [Fig.  4(a)].  Skaneateles  Lake 
had  the  lowest  minerogenic  particle  concentrations; 
these  were  particularly  lower  than  Owasco  Lake 
(next  lowest  A0  for  d  >  5  |xm. 

Seneca  River  N  levels  for  d  <  2  p,m  were  similar 
upstream  (B424)  and  downstream  (B409)  of  Cross 
Lake,  but  were  lower  for  larger  (>2  |xm)  minerogenic 
particles  at  the  downstream  site  [Fig.  4(b)].  Higher 
concentrations  were  observed  at  the  south  site  within 
Cross  Lake,  which  is  more  proximate  to  the  river 
inflow,  than  at  the  north  site.  A  substantial  decrease 
in  minerogenic  particle  concentrations  occurred  over 
the  river  reach  that  extends  from  just  downstream  of 
Cross  Lake  (B409)  to  a  position  18  km  further  down¬ 
stream  (B317). 

None  of  the  minerogenic  particle  populations  char¬ 
acterized  here  were  consistent  with  a  Junge  function. 
Application  of  this  function  resulted  in  overrepresen¬ 
tation  of  the  smaller  and  larger  particles  and  under¬ 
prediction  of  intermediate  size  particles  [an  example, 
the  south  site  of  Otisco  Lake,  is  illustrated  in  Fig.  4(c)]. 
The  values  of  j  ranged  from  2.5  to  3.2  (average  of  2.8), 
substantially  less  than  4,  a  value  commonly  invoked  in 
modeling  scattering  for  marine  systems  [5,35].  The 
2 C-B  performed  well  (r2  >  0.97  in  all  fittings)  in 
matching  the  observed  PSDs  in  the  fitted  size  range 
(i.e.,  >0.9  |xm),  and  in  most  cases  this  performance 
extended  to  smaller  sizes  [e.g.,  0.4-0. 5  |xm;  Fig.  4(c)]. 
With  the  exception  of  the  north  site  of  Otisco  Lake, 
values  of  jxfi  and  yB  (Table  5)  were  similar  to  the 
average  values  (2  and  0.226,  respectively)  reported  by 
Risovic  [4]  for  marine  particle  populations. 


Table  4.  Statistics  (Mean  ±  Standard  Deviation)  of  Minerogenic  Particle  Shapes  as  Described  by  ASP  (Aspect  Ratio)  for  Study  Sites 


Ca-Rich 

System  Diatom  Phacotus  Non -Phacotus  Ca-Agg  Clay  Quartz  Si- Rich 


Seneca  R. 


B424 

7.4 

± 

5.6 

— 

1.8 

± 

0.6 

2.4 

± 

2.2 

2.3  ± 

2.8 

1.6 

± 

1.0 

1.7  ±  0.5 

Cross  L.  North 

10.1 

± 

16.9 

1.1  ±  0.1 

1.7 

± 

0.8 

2.3 

± 

2.2 

2.1  ± 

1.2 

1.9 

± 

1.0 

1.8  ±  0.5 

Cross  L.  South 

6.6 

3.8 

1.1  ±  0.1 

1.7 

± 

0.6 

2.1 

± 

1.2 

2.3  ± 

2.0 

2.0 

± 

1.4 

1.7  ±  0.5 

B409 

7.2 

± 

5.6 

1.2" 

1.6 

± 

0.5 

2.1 

1.6 

2.2  ± 

1.8 

1.9 

± 

1.6 

1.7  ±  0.5 

B317 

5.7 

4.3 

— 

1.7 

± 

0.7 

2.0 

± 

1.4 

1.9  ± 

2.0 

1.5 

-t- 

0.4 

1.6  ±  0.4 

Lakes 

Owasco 

12.9 

± 

13.4 

1.2" 

2.3 

-f- 

2.6 

2.6 

-h 

5.0 

2.6  ± 

5.4 

2.7 

± 

2.6 

1.8  ±  0.5 

Skaneateles 

4.1 

± 

7.7 

— 

1.6 

0.8 

2.1 

± 

2.4 

2.2  ± 

3.0 

2.8 

± 

2.4 

1.7  ±  0.5 

Otisco  North 

3.4 

± 

0.7 

1.0" 

2.1 

± 

1.7 

2.3 

± 

1.7 

1.9  ± 

1.4 

1.4 

-+■ 

0.3 

1.5  ±  0.4 

Otisco  South 

6.0 

± 

7.7 

1.1  ±  0.0 

1.9 

± 

1.2 

2.2 

± 

1.8 

2.2  ± 

2.1 

1.6 

0.7 

1.6  ±  0.5 

Onondaga 

9.6 

+ 

8.4 

— 

1.5 

± 

0.4 

2.2 

-t- 

1.6 

2.4  ± 

2.4 

1.6 

± 

0.7 

1.6  ±  0.4 

"Less  than  three  observations. 
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Table  5.  2C -B  PSD  Fitting  Parameters 
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Fig.  4.  PSDs  for  minerogenic  particles:  (a)  five  lake  sites,  Owasco 
Lake,  Skaneateles  Lake,  north  and  south  sites  on  Otisco  Lake,  and 
Onondaga  Lake,  (b)  Seneca  River  system,  including  river  buoy  424, 
north  and  south  sites  on  Cross  Lake,  and  river  buoys  409  and  317, 
and  (c)  an  example  of  PSD  function  fits  [Junge  (j  =  2.97,  with  j  = 
4  added  for  reference)  and  2C-B;  OtS].  See  Table  1  for  site  abbre¬ 
viations. 


System" 

CB(1 

1  pm  l) 

Vb 

OtN 

5.012 

X 

1014 

5.25 

0.288 

Ow 

4.21 

X 

1014 

1.34 

0.172 

Sk 

1.58 

X 

1014 

2.64 

0.226 

OtS 

1.29 

X 

10 16 

2.82 

0.236 

On 

2.12 

X 

1015 

2.02 

0.199 

B424 

5.86 

X 

10'5 

1.62 

0.178 

CrN 

1.27 

X 

10'5 

1.54 

0.178 

CrS 

1.99 

X 

101S 

2.51 

0.217 

B409 

5.28 

X 

10,r> 

2.06 

0.202 

B317 

8.24 

X 

1014 

3.58 

0.262 

mean  ±  s.d/’ 

2.24  ±  0.72 

0.208  ±  0.03 

"See  Table  1  for  system  abbreviations. 

'’Mean  ±  standard  deviation,  OtN  not  included. 


someter  and  ac-s  (separate  casts)  approached  equiv¬ 
alence  [c(660)ac_s  =  1.09  X  c(660)r.Star  -  0.2]. 

Measured  values  of  6P(650)  ranged  from  0.47  m  1 
for  Skaneateles  Lake  to  7.9  m_1  at  B409  [Fig.  6(a)]. 
The  value  presented  for  the  south  site  of  Otisco  Lake, 
the  highest  of  the  study,  was  estimated  based  on 
laboratory  measurement  of  c(660)  (on  a  diluted  sam¬ 
ple)  according  to  fcp(650)  =  c(660)/1.045  [e.g.,  consis¬ 
tent  with  Babin  et  al.  [5]  for  6p(660)].  Four  of  the  sites 
had  6P(650)  values  of  ~4  m_1,  the  north  site  of  Otisco 
Lake,  Onondaga  Lake,  and  both  the  south  and  north 
sites  of  Cross  Lake.  Owasco  Lake,  Skaneateles  Lake, 
and  B317  had  6p(650)  values  <2  m'1. 

No  bhp(650)  observation  was  available  for  the  south 
site  of  Otisco  Lake  [Fig.  6(b)].  The  pattern  of  differ¬ 
ences  in  6/)p(650)  among  the  study  sites  was  strongly 
correlated  (r  =  0.96)  to  that  reported  for  6p(650). 
Values  of  6/,p(650)  ranged  from  0.0034  m  1  for  Ska¬ 
neateles  Lake  to  0.25  m  1  for  B424.  Despite  the 
strong  correlation  between  bp(6 50)  and  6,,p(650)  for 
this  population  of  sites,  substantial  differences  in  the 
backscattering  ratio  [66p(650)  =  6/,p(650)/6p(650)] 
were  observed  among  the  sites  [Fig.  6(c)].  The  ratio 


D.  Six  Hundred  Fifty  Nanometers  as  the  Reference 
Wavelength  for  b,  Bulk  Measurements  of  bp( 650)  and 
b6p(650),  and  Estimates  of  bm(650)  and  bb  m( 650) 

The  choice  of  \  =  650  nm  as  a  reference  wavelength 
(also  [39])  here  for  b  calculations  does  not  represent  a 
problem  for  comparing  the  results  from  this  study 
with  those  from  related  studies  that  have  adopted 
different  wavelengths.  For  example,  a  strong  (r2  = 
0.998),  nearly  1:1  (slope  =  1.04)  linear  relationship 
prevailed  between  6p(650)  and  the  spectrally  aver¬ 
aged  (400-730  nm)  scattering  coefficient  «6p);  Fig.  5). 
Similarly  strong  relationships  prevailed  between 
fcp(650)  and  other  selected  wavelengths,  including 
6P( 555)  (Fig.  5),  used  by  Babin  et  al.  [5]  as  a  reference 
wavelength.  The  values  of  c(660)  from  the  transmis- 


Fig.  5.  Relationships  between  (6P(\))  and  6P(650),  and  between 
6P(555)  and  6P(650),  for  study  sites,  with  equivalence  (1:1)  line  for 
reference. 
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Ow  Sk  UN  OlS  On  B424  CrN  CrS  B409  B317 


Fig.  6.  Comparisons  of  bulk  in  situ  measurements  of  scattering 
and  calculated  minerogenic  scattering  based  on  SAX  results:  (a) 
6P(650)  and  6m(650),  with  estimation  of  6„(650)  illustrated,  and  an 
example  of  b,n(650)  adjusted  for  the  lying  flat  effect  (inset),  (b) 
bhp{6 50)  and  bhm{6 50),  and  (c)  measured  and  minerogenic  6^(650). 
See  Table  1  for  site  abbreviations. 


was  <0.008  for  Skaneateles  and  Owasco  Lakes, 
—0.018  for  B424  and  —0.023  for  B317,  and  between 
0.0125  and  0.0145  for  the  remaining  sites. 


Estimated  values  of  6m(650),  bb >m(650),  and  the  cor¬ 
responding  minerogenic  66p(650)  are  presented  [Figs. 
6(a)-6(c)]  and  compared  to  the  overall  bulk  measure¬ 
ments.  Three  different  values  were  calculated  for 
6m(650)  and  6^(650)  for  each  site  [one  example  il¬ 
lustrated  in  inset  of  Fig.  6(a)]  that  correspond  to  Mie 
theory  calculations  based  on  SAX  results:  (1)  without 
adjustments  for  the  particle  “lying  flat”  effect,  (2) 
with  a  reasonable  lower  bound  representation  of  the 
effect  (ASP  >  2.5),  and  (3)  with  a  reasonable  upper 
bound  representation  of  the  effect  (ASP  >  1.5).  The 
adjustments  resulted  in  decreases  in  the  estimated 
values  of  fe,„(650)  [inset  of  Fig.  6(a)]  and  bh  m(650).  The 
minimum  decreases,  for  the  lower  bound  adjust¬ 
ments,  was  about  4%  for  6m(650)  and  3 %  for  bh  m(650); 
the  averages  were  about  5%  for  both  features  of  light 
scattering.  The  maximum  decreases,  for  the  upper 
bound  adjustments,  were  —15%;  the  averages  were 
12.6  and  13.5%  for  6,„(650)  and  bhm(6 50),  respec¬ 
tively.  The  single  best  estimate  for  each  site  is  pre¬ 
sented  here  [Figs.  6(a)  and  6(b)]  as  the  average  of  the 
values  for  upper  and  lower  bound  adjustments.  These 
values  have  been  used  in  evaluations  of  relationships 
with  measurements  and  in  estimation  of  6(/650)  val¬ 
ues  [Fig.  6(a)]. 

Differences  in  ( Qbm )  (Table  6)  were  modest  between 
the  sites  and  among  the  minerogenic  classes  for  in¬ 
dividual  sites  (Table  6);  values  of  ( Qb  m )  ranged  from 
2.1  to  2.66  for  the  sites.  The  clay  mineral  values  for 
( Qbm )  were  most  often  the  lowest  among  the  classes; 
the  average  for  this  class  for  the  10  sites  was  2.23. 
The  average  for  quartz  (2.43),  a  relatively  minor  com¬ 
ponent  of  bmi  was  9%  higher.  The  overall  average 
( Qbb,m )  for  the  study  sites  was  0.061,  with  a  coefficient 
of  variation  of  7.6%. 

Comparison  of  calculated  6,„(650)  and  bbm(650) 
with  bulk  in  situ  measurements  showed  that  the  min¬ 
erogenic  component  was  a  substantial  contributor  to 
total  scattering  and  backscattering.  The  minerogenic 
component  of  bp(6 50)  ranged  from  —20%  at  Owasco 
Lake  to  94%  at  B424.  This  component  represented 
from  30%  to  50%  of  6P(650)  at  six  of  the  sites;  the 
average  for  all  sites  was  42%.  These  contributions 


Table  6.  Mean  Scattering  Efficiencies,  (Qb  m(650)>,"  of  Minerogenic  Particle  Types  for  Study  Sites 


System 

Clay 

Ca-Rich 

Non -Phacotus 

Quartz 

Si-Rich 

Ca-Agg 

Miscellaneous 

Total 

Seneca  R. 

B4‘24 

2.15 

2.23 

2.4 

2.28 

2.14 

2.48 

2.18 

Cross  L.  North 

2.29 

2.19 

2.44 

2.15 

2.23 

2.17 

2.22 

Cross  L.  South 

2.24 

2.12 

2.28 

2.12 

2.25 

2.37 

2.22 

B409 

2.19 

2.43 

2.4 

2.14 

2.27 

2.32 

2.23 

B317 

2.26 

2.33 

2.53 

2.36 

2.43 

2.48 

2.33 

Lakes 

Owasco 

2.18 

2.25 

2.4 

2.25 

2.23 

2.24 

2.23 

Skaneateles 

2.26 

2.46 

2.43 

2.28 

2.25 

2.43 

2.3 

Otisco  North 

2.15 

2.28 

2.58 

2.53 

2.12 

2.1 

2.2 

Otisco  South 

2.3 

2.66 

2.41 

2.18 

2.38 

2.13 

2.3 

Onondaga 

2.3 

2.2 

2.42 

2.42 

2.18 

2.39 

2.28 

nFor  case  of  lying-flat  adjustments  for  particles  with  ASP  >  1.5. 
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were  shifted  higher  for  backscattering.  Estimates  of 
6ft  m(650)  represented  between  50%  and  80%  of  the 
bbp(6 50)  measurements  at  eight  of  nine  sites;  the  av¬ 
erage  for  all  sites  was  74%.  The  bb m(650)  estimate  for 
B424  exceeded  the  6>,p(650)  observation  by  30%.  Es¬ 
timates  of  the  bbftn( 650)  varied  from  —0.25  for  the 
river  sites  to  —0.32  for  Owasco  Lake  and  the  north 
site  on  Otisco  Lake  [Fig.  6(c)].  These  ratio  values  in 
all  cases  exceeded  those  based  on  measurements,  by 
factors  of  more  than  two,  except  for  the  river  sites. 

E.  Relationships  between  Measures,  Estimates,  and 
Indicators  of  Scattering 

TSS  was  a  strong  predictor  of  the  differences  of  both 
6P(650)  and  bbp(6 50)  [(Figs.  7(a)  and  7(b)]  among  the 
study  sites,  explaining  97%  (p  <  0.001)  and  94% 


( p  <  0.001)  of  the  differences,  respectively,  according 
to  linear  least-squares  regression  analysis.  The  mean 
values  of  fep(650)/TSS  and  6/,p(650)/TSS  ratios  were 
0.86  and  0.012  m2g_1,  respectively.  PAV,„  was  a 
strong  predictor  of  the  observed  differences  in  6p(650) 
and  6/,p(650)  [(Figs.  7(c)  and  7(d)],  explaining  96% 
(p  <  0.001)  of  the  differences  in  both  cases,  a  result 
consistent  with  the  performance  of  FSScorr  in  predict¬ 
ing  c(660)  [Eq.  (9)]. 

FSScorr  was  nearly  as  strong  a  predictor  of  the  es¬ 
timated  values  of  6m(650)  and  66  „,(650)  [Figs.  7(e)  and 
7(f)]  as  TSS  was  for  the  overall  bulk  scattering  met¬ 
rics,  accounting  for  95%  (p  <  0.001)  of  the  variations 
in  both  estimates.  The  mean  values  of  the 
6m(650)/FSScorr  and  bh  m(650)/FSScorr  ratios  were  0.62 
and  0.017  nr  g  l,  respectively.  Significant  relation- 
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Fig.  7.  Evaluation  of  relationships  between  measures,  estimates,  and  indicators  of  scattering  for  the  study  sites:  (a)  6^(650)  versus  TSS, 
(b)  6^(650)  versus  TSS,  (c)  6P(650)  versus  PAV,„  (B424  as  outlier  in  regression  analysis),  (d)  6^(650)  versus  PAV,„  (£424  as  outlier  in 
regression  analysis),  (e)  6„,(650)  versus  FSScorr,  (f)  66>„,(650)  versus  FSScorr,  (g)  bn ,(6501/6^(650)  versus  FSScon/TSS,  (h)  66t>n(650)/6fc/,(650) 
versus  FSScorl./TSS,  (i)  66/}(650)  versus  FSScorr/TSS,  (j)  bJ650)  versus  VSScolT,  (k)  6„(650)  versus  PAV„,  and  (1)  6o(650)  versus  6oChi(650).  See 
Table  1  for  site  abbreviations. 
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ships  were  observed  between  the  FSScorr/TSS  ratio 
and  the  6m(650)/6p(650)  ratio  [Fig.  7(g)]  as  well  as  the 
66m(650)/fc/)p(650)  ratio  [Fig.  7(h)].  However,  the  ob¬ 
servations  for  B424  had  a  disproportionate  effect  in 
both  cases;  the  relationships  were  not  significant  if 
the  values  for  that  site  were  omitted,  and  in  the  case 
of  bhmy  the  estimate  exceeded  the  measured  [Fig. 
7(h)].  Some  positive  dependence  between  6^(650) 
and  FSScorr/TSS  was  observed  across  the  study  sites 
[Fig.  7(i)];  however,  even  with  B317  as  an  outlier, 
the  relationship  was  not  significant  at  a  =  0.05 
(p  =  0.07). 

A  strong  linear  relationship  (r2  =  0.89,  p  <  0.001), 
with  a  relatively  small  y  intercept,  prevailed  between 
feo(650)  [=6P(650)  -  6m(650)]  and  VSScorr  across  the 
study  sites  [Fig.  7(j)].  The  mean  fto(650)/VSScorr  ratio 
was  1.08  m2  g1  for  the  study  sites.  The  relationship 
observed  between  6o(650)  and  the  PAV  associated 
with  organic  particles  [PAVf„  Fig.  7(k);  estimated  as 
the  sum  of  Phacotus  and  the  Organics  and  Diatoms 
classes,  without  lying  fiat  adjustments]  was  nearly  as 
strong.  The  best  fit  relationship  between  calculated 
6o(650)  and  the  independently  estimated  6oChl(650) 
approached  1:1,  though  it  was  less  strong  (r2  = 
0.53,  p  =  0.02)  [Fig.  7(1)]. 

F.  Size  Dependencies  of  Scattering  and  Sensitivity  to  n 
Value 

Calculated  (Mie  theory)  cumulative  size  distributions 
of  6/n(650)  and  bh  m(6 50),  based  on  SAX  observations, 
are  presented  to  depict  the  contributions  of  various 
particle  sizes  for  a  Seneca  River  (B424)  site  [Figs.  8(a) 
and  8(b)].  These  patterns  are  generally  representa¬ 
tive  of  those  obtained  for  all  the  study  sites.  Particles 
within  the  size  range  2-10  jxm  contributed  the  most 
to  both  6;n(650)  and  bhm(6 50).  Larger  particles  were 
generally  unimportant  due  to  their  scarceness,  as 
depicted  by  the  “leveling  off’  of  plots  for  d  >  —10  |xm. 
Submicron  particles  made  <2%  contributions  to  scat¬ 
tering.  The  size  contribution  patterns  were  generally 
similar  for  6m(650)  and  fe6m(650).  The  size  that  cor¬ 
responds  to  50%  of  the  total  value  (i.e.,  particles  less 
than  or  equal  to  that  size  are  responsible  for  50%  of 
the  scattering;  d50)  is  a  useful  attribute  of  the  pat¬ 
terns.  The  db0  values  for  both  6,n(650)  and  56m(650) 
for  the  river  site  were  ~5.2  |xm. 

The  size  dependency  patterns  of  scattering  were 
also  simulated  with  uniform  m  (1.17-0.005D  for 
three  cases  of  PSDs:  observed,  best  fit  Junge,  and 
Junge  with  j  =  4  (curves  in  Fig.  8),  and  the  simula¬ 
tions  for  the  Junge  cases  were  extended  down  to  a 
size  of  0. 1  |jtm.  The  simulation  performed  for  the  mea¬ 
sured  PSDs  (with  uniform  m)  tracked  the  observed 
patterns  closely.  Consistent  with  the  shortcomings  of 
the  best  fit  Junge  function  in  fitting  the  observed 
PSDs  [e.g.,  Fig.  4(c)],  this  function  overrepresented 
the  scattering  effects  of  submicrometer  particles  and 
predicted  continuing  strong  increases  in  scattering 
for  d  >  10  |xm,  and  generally  deviated  from  observa¬ 
tions  for  intermediate  sizes  (Fig.  8).  The  deviations  of 
the  shapes  of  the  patterns  for  j  =  4  were  even  greater 


Fig.  8.  Size  dependencies  of  scattering  coefficients  as  cumulative 
plots,  based  on  observations  (symbols)  from  Seneca  River  at  buoy 
424,  and  simulations  from  three  cases  of  PSDs  (observed,  and 
Junge  with  best  fitted  j  and  invoked  j  =  4;  curves)  with  uniform 
m  =  1.17-0.0005*:  (a)  bJ6 50),  and  (b)  6,,„,(6 50). 

(for  the  case  of  maintaining  the  same  C  value;  Fig.  8). 
In  contrast,  the  2C-B  fits  supported  cumulative  scat¬ 
tering  patterns  (not  shown)  that  tracked  those  based 
on  observations  closely,  consistent  with  their  good  fits 
of  the  observed  PSDs  in  the  size  range  that  was  im¬ 
portant  to  scattering. 

Contrasting  patterns  emerged  for  the  dependen¬ 
cies  of  6m(650)  and  6/,m(650)  on  n,  as  illustrated  for 
the  Owasco  Lake  sample  [Figs.  9(a)  and  9(b)].  Values 
are  presented  relative  to  the  predictions  for  n  =  1.25 
to  support  comparisons  for  the  observed  PSDs  and 
Junge  function  representations.  The  predicted  pat¬ 
terns  for  this  site  were  qualitatively  representative  of 
the  others  included  in  this  study.  Values  of  6„,(650) 
were  predicted  to  be  relatively  invariant  over  the 
evaluated  range  of  n  (1.05-1.25)  for  the  observed 
PSDs  [Fig.  9(a)],  whereas  progressive  increases  in 
bb m(650)  were  predicted  throughout  the  range  for  the 
observed  PSDs  [Fig.  9(b)].  Junge  function  scenarios 
for  the  PSDs  deviated  from  the  observed,  demonstrat¬ 
ing  a  greater  positive  dependence  on  n,  with  an  even 
stronger  effect  for  j  =  4  compared  to  the  best  fit  value 
of  2.68. 

4.  Discussion 

A.  Importance  and  Partitioning  of  the  Minerogenic 
Particle  Assemblage 

The  central  role  minerogenic  particles  played  in  reg¬ 
ulating  differences  in  light-scattering  levels  between 
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Fig.  9.  Scattering  as  functions  of  n  (with  n'  =  0.0005),  for  ob¬ 
served  minerogenic  PSD  at  Owasco  Lake,  for  Junge  PSDs  with 
fitted  j  values  of  2.68  and  commonly  invoked  4:  (a)  6„,(650)  as 
percentages  of  the  values  calculated  for  n  =  1.25  (0.41,  0.45,  and 
0.22  m  1  for  the  three  cases  of  PSDs,  respectively),  and  (b)  bhm(6 50) 
as  percentages  of  the  values  calculated  for  n  =  1.25  (0.016,  0.016, 
and  0.0092  m  1  for  the  three  cases  of  PSDs,  respectively). 


the  study  systems  was  manifested  in  empirical  anal¬ 
yses  with  common  metrics  of  water  quality  [Eq.  (9)], 
as  well  as  by  the  dependence  of  both  6P(650)  and 
bbJ6 50)  on  PAVm  [Figs.  7(c)  and  7(d)].  Increasingly, 
lakes  and  reservoirs  have  been  identified  where  min¬ 
erogenic  particles  are  important  in  regulating  scat¬ 
tering  and  apparent  optical  properties  (AOPs).  For 
example,  Peng  and  Effler  [27]  reported  that  PAVm 
was  a  stronger  (r2  =  0.82)  predictor  than  Chi  (r2 
=  0.56)  of  Secchi  depth  differences  between  the 
eleven  Finger  Lakes  of  New  York  (see  also[29,40,41]). 
Moreover,  distinct  temporal  structure  in  optical  char¬ 
acteristics  in  freshwater  systems  can  be  imparted  by 
dynamics  in  terrigenous  inputs  [15]  and  autochtho¬ 
nous  precipitation  of  minerogenic  particles  [29,42]. 
The  origin  of  optically  important  particles  (allochtho¬ 
nous  versus  autochthonous)  is  essential  information 
for  managers  concerned  with  identifying  appropriate 
targets  and  rehabilitation  approaches  for  improved 
optical  esthetics,  and  for  proper  interpretation  of  bulk 
optical  measurements. 

Despite  the  inherent  limitations  of  a  single  sam¬ 
pling,  the  PAVm  results  are  generally  consistent 
with  known  characteristics  of  the  study  systems. 
All  of  the  study  lakes  are  known  to  be  oversaturated 
with  respect  to  calcite,  and  to  experience  precipita¬ 
tion  of  CaC03  in  the  upper  waters  during  summer 
[29,42,43].  The  spatial  PAY  patterns  for  clay  mineral 


particles  and  diatoms  from  upstream,  within,  and 
just  downstream  of  Cross  Lake  reflect  the  incomplete 
mixing  of  the  river  inflow  within  this  lake  [16]  and 
the  functioning  of  the  lake  as  a  source  of  phytoplank¬ 
ton  [e.g.,  diatoms;  see  also  Fig.  2(c)]  to  downstream 
portions  of  the  river  [17].  The  decrease  in  PAVm  (Ta¬ 
ble  3)  and  FSScorr  [Fig.  2(b)]  between  B409  and  B317 
reflects  the  effects  of  the  nonselective  filter  feeding  by 
zebra  mussels  [44];  this  sink  for  inorganic  (as  well  as 
organic)  particles  has  contributed  importantly  to  the 
increased  clarity  observed  over  this  reach  since  the 
invasion  [45]  by  this  exotic  bivalve. 

This  paper  has  embraced  the  reductionist  approach 
described  by  Stramski  and  co-workers  [7,8]  by  parti¬ 
tioning  the  minerogenic  fraction  of  scattering  into 
multiple  components.  SAX  could  support  partitioning 
of  the  minerogenic  assemblage  into  more  particle 
classes  than  the  number  adopted  here.  For  example, 
Yin  and  Johnson  [46]  used  a  particle  classification 
scheme  that  delineated  18  minerogenic  classes,  in¬ 
cluding  five  types  of  clay  minerals,  in  applying  SAX 
for  a  particle-class  budget  analysis  of  Onondaga  Lake 
sediments.  Such  schemes  use  tighter  boundaries  for 
the  elemental  stoichiometries  of  the  particle  classes, 
and  implicitly  are  accompanied  by  greater  uncer¬ 
tainty  with  regards  to  the  proper  classification  of  in¬ 
dividual  particles.  While  such  an  approach  may  be 
appropriate  for  specific  issue/ecosystem  combina¬ 
tions,  the  simpler  scheme  adopted  here  (Table  2) 
served  to  resolve  the  fundamentally  different  origins 
of  the  important  components  of  the  minerogenic  par¬ 
ticle  assemblages  of  these  study  systems.  Further 
compression  of  the  terrigenous  particle  types  into 
fewer  classes  is  not  desirable,  even  though  the  n  val¬ 
ues  are  similar,  as  differences  in  n '  may  have  impor¬ 
tant  implications  with  respect  to  remote-sensing 
reflectance  spectral  signatures  [35]. 

B.  Implementation  of  the  Forward  Method, 

Partitioning  b(\) 

Advancement  of  the  reductionist  approach  to  par¬ 
tition  scattering  requires  implementation  of  the 
forward  method,  the  calculation  of  scattering  at¬ 
tributable  to  noteworthy  particle  classes  based  on 
particle  information.  Limitations  in  our  implemen¬ 
tation  of  the  forward  method  have  two  sources:  (1) 
uncertainties  in  the  information  supplied  to  the  calcu¬ 
lation  framework,  and  (2)  application  of  the  Mie  theory 
for  the  documented  particle  populations.  The  capabil¬ 
ity  of  SAX  to  characterize  individual  minerogenic  par¬ 
ticles  both  compositionally  and  morphometrically 
supports  the  advancement  of  the  forward  method. 
There  are  modest  imperfections  in  both  regards  rela¬ 
tive  to  supporting  forward  calculations.  For  example, 
there  are  only  minor  uncertainties  in  the  partitioning 
of  minerogenic  particles  for  the  simple  classification 
strategy  adopted  because  of  the  conspicuous  elemental 
differences  (Table  2).  Given  the  elemental  composition 
information  from  SAX,  the  greater  uncertainty  lies  in 
the  specification  of  n  values  based  on  the  literature. 
Noteworthy  variations  in  n  are  associated  with  the 
various  clay  minerals  and  the  source  of  the  value  for 
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CaC03  (calcite)  (1.24  by  Twardowski  et  al.  [47]  and 
1.173  by  Wozniak  and  Stramski  [35]).  Estimates  of 
6m(650)  were  insensitive  to  minerogenic  n  values  for 
the  PSDs  encountered  at  the  study  sites  [Fig.  9(a)]. 
However,  estimates  of  bh >m(650)  were  sensitive  to  this 
source  of  uncertainty  [Fig.  9(b)],  indicating  the  im¬ 
portance  of  the  composition  information  from  SAX  in 
guiding  specification  of  n  in  calculations  of  backscat- 
tering.  SAX  provides  a  robust  representation  of  par¬ 
ticle  PA.  The  shortcoming  of  lying  flat  on  the  filter  in 
SAX  is  modest  compared  to  that  associated  with 
Coulter  counter  data  for  irregularly  shaped  particles. 
Such  counters  infer  particle  sizes  from  measured 
electrical  resistance  of  particles  based  on  their  vol¬ 
umes,  leading  to  systematic  underestimation  of 
particle  sizes  (calculated  from  volume  equivalent 
spheres)  for  nonspherical  particles.  The  presented 
bounds  for  the  lying  flat  effect  are  believed  to  reason¬ 
ably  bracket  the  effect. 

The  Mie  theory  stipulations  of  particle  sphericity 
and  homogeneity  were  not  met  for  the  minerogenic 
particles  of  this  study,  as  reflected  by  the  ASP  values 
(Table  4)  and  substantial  concentrations  of  aggre¬ 
gates  of  minerogenic  classes  (e.g.,  Ca-agg;  Table  3). 
These  deviations  need  to  be  considered  in  context. 
First,  the  heterogeneity  associated  with  partial  coat¬ 
ing  of  particles  with  CaC03  is  a  minor  effect  given  the 
similarity  of  n  values  for  the  various  classes  (Table  2) 
and  the  insensitivity  of  6m(650)  estimates  to  varia¬ 
tions  in  n  (Fig.  9).  Second,  it  has  long  been  recognized 
that  most  oceanic  (and  presumably  freshwater)  par¬ 
ticles  are  not  spherical  [14].  Deviations  from  spheric¬ 
ity  can  cause  substantial  shifts  in  Pp(0),  though  these 
effects  influence  light  scattered  in  the  backward  di¬ 
rection  more  than  in  the  forward  (and  primary)  di¬ 
rection  [6,48,49].  Both  positive  and  negative  effects  of 
deviations  from  sphericity  on  bhp{\)  have  been  re¬ 
ported  [50,51].  The  effect  of  particle  shape  for  ran¬ 
domly  oriented  particles  on  ( Qb )  (and  consequently  bp) 
is  apparently  small  by  comparison  [14,48,52].  Thus 
we  have  relatively  more  confidence  in  our  Mie  theory 
estimates  of  6m(650)  than  those  of  bbm(6 50).  Mie  the¬ 
ory  calculations  represent  the  only  viable,  albeit  im¬ 
perfect,  vehicle  to  estimate  both  bm(\)  and  66m(X) 
based  on  individual  particle  information  for  thou¬ 
sands  of  particles  that  differ  in  size  and  composition 
[53].  More  complex  modeling,  if  feasible  for  such  het¬ 
erogeneous  particle  populations,  would  still  be  ac¬ 
companied  by  substantial  uncertainty  that  is  difficult 
to  assess  [53].  Minerogenic  particles  such  as  those 
encountered  in  our  study  (e.g.,  Fig.  3)  represent  less 
complex  problems  in  light  scattering  compared  to 
many  microorganisms  with  respect  to  geometry  [51] 
and  the  potential  for  layered  effects  within  microor¬ 
ganisms  [8],  though  variation  in  the  degree  of  hydra¬ 
tion  [20]  may  be  an  issue  (albeit  modest)  for  clay 
minerals. 

Despite  the  deviations  of  the  particle  populations 
from  the  simplifying  assumptions  of  Mie  theory,  the 
evidence  is  strong  that  SAX  can  support  the  parti¬ 
tioning  of  the  minerogenic  component  of  scattering  in 
these  fresh  waters  according  to  the  relative  contribu¬ 


tions  of  important  particle  classes.  This  partitioning 
can  reasonably  be  based  directly  on  PAVm  results 
delineated  according  to  designated  classes  (Table  3). 
This  position  is  supported  by  the  general  similarities 
of  (Qb,m)  (Table  6;  embedded  effects  of  m  and  PSD) 
and  ASP  representations  of  shape  for  the  various 
classes  (Table  4).  Accordingly,  percent  contributions 
of  the  various  minerogenic  particle  classes  to  PAVm 
correspond  to  contributions  to  bm\  e.g.,  clay  minerals 
represented  —70%  of  6m(650)  at  the  upstream  river 
site  (B424)  and  calcite  (non -Phacotus  Ca-rich)  was 
responsible  for  —35%  of  6m(650)  in  Owasco  Lake 
(Table  3). 

Demonstration  of  closure  of  the  SAX-based  esti¬ 
mates  with  independent  measures  of  minerogenic 
scattering  is  difficult  for  natural  waters  because 
particle  populations  are  composed  of  mixtures  of 
minerogenic  and  organic  particles.  The  closure 
between  SAX-based  estimates  of  bm(66 0)  and  mea¬ 
surements  of  c(660)  demonstrated  for  Schoharie  Res¬ 
ervoir,  where  clay  particles  dominate  scattering  [15], 
is  supportive  of  the  minerogenic  scattering  estimates 
presented  here.  None  of  the  sites  of  this  study  dem¬ 
onstrated  the  extent  of  dominance  by  minerogenic 
particles  [e.g.,  Fig.  2(b)]  reported  for  Schoharie  Res¬ 
ervoir.  However,  the  significant,  nearly  1:1  relation¬ 
ship  between  the  calculated  6o(650)  values  and 
independent  estimates  of  b()  Chi(650)  suggests  a  degree 
of  closure.  Much  of  the  observed  scatter  in  the 
relationship  may  reflect  the  limitations  of  Chi  as  a 
surrogate  metric  of  the  overall  organic  particle  pop¬ 
ulation.  Moreover,  the  array  of  measured  and  esti¬ 
mated  conditions  included  in  this  study  offered 
multiple  tests  of  consistency  that  generally  support 
the  SAX-based  estimates  of  minerogenic  scattering. 
This  is  supported  by  the  strong  relationships  between 
the  gravimetric  measure  of  minerogenic  suspended 
solids  (FSScorr)  and  fe„,(650)  [Fig.  7(e)],  and  between 
FSScorr  and  bb/n(650)  [Fig.  7(f)].  Further,  the  ratios  of 
6m(650)  to  6p(650),  and  6^(650),  were  positively  re¬ 
lated  to  the  FSScorr/TSS  [Figs.  7(g)  and  7(i)].  Addi¬ 
tional  support  for  the  6m(650)  predictions  is  provided 
by  the  consistency  of  the  ba(6 50)  estimates  with  indi¬ 
cators  of  the  organic  particle  population,  including 
the  strong  relationships  between  6o(650)  and  VSScorr 
[Fig.  7(j)],  and  between  feo(650)  and  PAV„  [Fig.  7(k)j. 
Overall,  the  empirical  indicators  of  the  performance 
of  SAX  in  supporting  the  partitioning  of  bp(\)  into 
bm(\)  and  60(\)  are  promising.  We  have  not  pursued 
similar  partitioning  for  bhp(\)  here  because  of  the 
greater  uncertainty  associated  both  with  the  Mie  es¬ 
timates  of  bh  m(\)  and  empirical  relationships  to  pre¬ 
dict  bbo(\)  based  on  Chi  [38]. 

The  extent  of  consistency  of  the  scattering  estimates 
depends  importantly  on  the  representativeness  of  the 
samples  relative  to  the  field  instrumentation  measure¬ 
ments.  The  general  consistency  of  the  6p(650)/TTS 
ratios  observed  here  with  the  marine  literature  [5] 
offers  only  coarse  support  for  the  representative¬ 
ness  of  the  samples.  Space-time  differences  in  the 
waters  sampled  and  analyzed  in  the  laboratory  ver- 
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sus  measured  in  situ ,  associated  with  patchy  distri¬ 
butions  of  particles  (e.g.,  Babin  et  al.  [5]),  likely 
contributed  to  the  observed  scatter  in  the  relation¬ 
ships  (Fig.  7)  and  the  inconsistent  case  of  66m(650) 
>  bhp(650)  at  B424  [Fig.  6(b)].  The  modest  variations 
in  the  relationship  between  the  ac-s  and  C-Star  mea¬ 
surements  of  c(660)  among  the  sites  suggests  that 
such  monitoring  limitations  likely  contributed  to  im¬ 
perfect  relationships. 

C.  Minerogenic  Particle  Size  Distributions 

Presentation  of  detailed  minerogenic  PSDs  that  in¬ 
cluded  submicrometer  particles  (e.g.,  Fig.  4)  has  been 
rare.  Lambert  et  al.  [54]  used  scanning  electron  mi¬ 
croscopy  integrated  with  an  electron  microprobe  tech¬ 
nique  to  characterize  the  PSDs  of  minerogenic 
particles  >0.2  p.m  in  15  samples  from  seven  ocean 
sites.  They  reported  a  leveling-off  of  particle  concen¬ 
trations  for  sizes  <0.7  p.m,  deviating  strongly  from  a 
Junge  function,  and  represented  the  PSDs  by  lognor¬ 
mal  distributions  instead.  The  PSDs  reported  here 
(Fig.  4)  had  similar  patterns  to  those  reported  re¬ 
cently  [15]  for  Schoharie  Reservoir,  based  on  SAX 
measurements;  e.g.,  conspicuous  decreases  in  F(d)  for 
sizes  smaller  than  the  frequency  peak  size  (exception 
for  north  site  of  Otisco  Lake)  within  the  submicrome¬ 
ter  range.  We  are  unaware  of  any  direct  measure¬ 
ments  that  included  submicrometer  observations  for 
either  marine  or  fresh  waters  to  support  the  Junge 
function  as  representative  of  minerogenic  PSDs. 
Coulter  counters  cannot  address  this  issue  because 
the  lower  size  threshold  is  too  high,  and  these  instru¬ 
ments  fail  to  differentiate  minerogenic  from  organic 
particles. 

Characterization  of  PSDs  has  received  extensive 
attention  for  marine  systems  [2,55].  A  number  of 
functional  representations  of  PSDs  have  been  used, 
including  the  hyperbolic  (or  Junge)  model,  the  lognor¬ 
mal  distribution  [54],  segmented  lognormal  functions 
[2],  the  TCM,  and  segmented  hyperbolic  distributions 
[55].  The  Junge  function  was  supported  in  earlier  stud¬ 
ies  that  used  electronic  particle  counters  (e.g.,  [56]), 
and  has  been  widely  adopted  in  modeling  [1,5,7]  and 
inversion  analyses  [47,57].  An  exception  was  the  inclu¬ 
sion  of  the  TCM,  in  addition  to  a  Junge  function,  in  a 
modeling  analysis  by  Twardowski  et  al.  [47].  Further 
evaluation  is  needed  to  test  the  applicability  of  2C -B  as 
well  as  other  functions  in  representing  PSDs  in  other 
fresh  waters.  Other  models  may  also  fit  the  observa¬ 
tions,  particularly  the  lognormal  function.  The  gener¬ 
ality  of  such  models  is  important  to  support  inversion 
algorithms. 

Babin  et  al.  [5]  predicted  submicrometer  mineral 
particles  (n  =  1.18)  contribute  substantially  to  bm  in 
marine  waters,  with  a  value  of  d50  ~  1  fxm  for  j  =  4 
based  on  Mie  theory  calculations  (our  simulations  of 
the  size  dependency  patterns  of  scattering  adopting 
j  =  4  resulted  in  a  similar  conclusion;  dotted  curves  in 
Fig.  8).  Accordingly,  those  investigators  suggested  that 
the  determination  of  the  bp/ SPM  ratio  may  be  sensi¬ 
tive  to  the  porosity  of  the  filter  used  in  their  SPM 
analyses  (glass  fiber  GF/F  filters,  0.7  |xm  pore  size). 


Our  observations  for  these  freshwater  systems  are 
strikingly  different  in  that  regard.  Submicrometer 
minerogenic  particles  did  not  make  noteworthy  contri¬ 
butions  to  scattering  according  to  Mie  theory  calcula¬ 
tions  based  on  SAX  data  (Fig.  8);  e.g.,  d50  ~  4  p.m. 
Moreover,  the  effective  pore  size(s)  of  glass  fiber  filters 
have  been  reported  to  be  less  than  the  nominal  sizes 
[58,59]  (e.g.,  effective  size  of  our  filters  <1.5  pm;  a 
better  choice  would  have  been  GF/F  with  a  nominal 
pore  size  of  0.7  pm).  Under  these  conditions,  the 
effective  pore  size  of  the  filters  used  in  the  TSS  an¬ 
alyses  does  not  represent  a  noteworthy  source  of 
uncertainty  for  the  fem(650)/FSScorr  ratio.  We  cannot 
eliminate  the  possibility  that  large  quantities  of  sub¬ 
micrometer  organic  particles  were  present  and  passed 
through  the  filters  and  thereby  caused  the  6o(650)/ 
VSS C01T  and  6/;(650)/TSS  ratios  to  be  false  high.  How¬ 
ever,  review  of  the  available  particle  size  information 
for  the  Organics  class  for  the  study  sites  suggests  sub¬ 
micron  organic  particles  were  not  important. 

The  limitations  in  the  best  fit  J  function  in  match¬ 
ing  the  observed  PSDs  of  this  study  resulted  in  note¬ 
worthy  systematic  shortcomings  in  representation  of 
the  dependence  of  6m(650)  and  6/)  m(650)  on  particle 
size  (Fig.  8).  For  example,  though  the  predicted 
6m(650)  approximately  matched  that  based  on  obser¬ 
vations  through  a  size  range  extending  to  —20  pm, 
the  size  trajectory  of  the  relationships  deviated  sub¬ 
stantially,  and  unrealistic  strong  increases  beyond 
the  sizes  of  SAX  observations  (i.e.,  >20  pm)  were 
predicted  [Fig.  8(a)].  Similar  shortcomings  of  the  J 
function  in  representing  the  size  dependence  of  bm 
have  been  reported  by  Stavn  and  Keen  [60].  More¬ 
over,  invoking  a  J  functional  representation  of  the 
minerogenic  PSDs  of  this  study  caused  the  depen¬ 
dence  of  calculations  of  b„ ,(650)  and  bhm(6 50)  on  n 
(i.e.,  particle  composition)  to  be  overstated  (Fig.  9). 

D.  Backscattering 

Great  uncertainty  persists  concerning  the  origins  of 
backscattering  in  both  marine  [8,38,39]  and  freshwa¬ 
ter  [9]  systems,  a  key  issue  in  remote  sensing  [50]. 
Stramski  et  al.  [8]  described  the  dominant  fraction  of 
backscattering  that  was  not  accounted  for  by  micro¬ 
organisms  as  the  “missing  backscattering  enigma.” 
Submicrometer  nonliving  (particularly  inorganic)  par¬ 
ticles  were  identified  as  the  likely  responsible  compo¬ 
nent  in  the  open  ocean  under  nonbloom  conditions 
because  these  particle  sizes  make  relatively  greater 
contributions  to  bhp{\)  than  to  bp(\)  [1,7].  Similarly, 
modeling  analyses  by  Wozniak  and  Stramski  [35]  de¬ 
picted  the  potential  importance  of  submicrometer 
minerogenic  particles  in  influencing  ocean  reflec¬ 
tance.  Modeling  analyses  by  Stramski  et  al.  [7],  based 
on  speculative  concentrations  and  PSDs  (j  =  4),  sug¬ 
gested  that  submicrometer  minerogenic  particles 
were  important  for  bp{\)  and  dominated  bhp(\)  in 
ocean  waters.  The  argument  that  backscattering  is 
largely  regulated  by  submicrometer  particles  has 
been  influenced  greatly  by  the  lack  of  direct  measure¬ 
ments  of  these  colloidal  particles,  which  has  necessi¬ 
tated  extrapolation  of  fitted  or  assumed  functions 
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into  this  size  range.  Risovic  [4]  demonstrated  that 
bbp(\)  depended  strongly  on  the  trajectory  of  the  ex¬ 
trapolation  and  thus  the  specified  relationship  to  de¬ 
scribe  the  PSD.  Based  on  comparative  fits  of  the 
Junge  and  the  TCM  to  observed  ocean  PSDs  (d 
>  0.8  |xm),  Risovic  [4]  demonstrated  that  application 
of  the  Junge  function  resulted  in  substantial  overpre¬ 
diction  of  the  bbp(\),  as  a  result  of  overestimation  of 
concentrations  of  submicrometer  particles. 

Several  features  of  the  results  established  that 
submicrometer  minerogenic  particles  were  not  im¬ 
portant  contributors  to  bbp(6 50):  (1)  the  observed 
PSDs  were  fitted  well  with  the  2 C-B  function,  (2)  the 
magnitude  of  bbm(6 50)  relative  to  the  bhp(6 50)  obser¬ 
vations  [Fig.  6(b)],  (3)  the  size  trajectories  of  bbm( 650) 
[Figs.  8(b)  and  8(d)],  and  (4)  the  lack  of  dependence 
of  ASP  on  d.  The  case  of  invoking./  =  4  for  the  PSDs 
of  this  study  demonstrates  the  preferential  effect  on 
bbni(65 0)  versus  6m(650)  for  dependencies  on  d  (Fig. 
8).  Accordingly,  —60%  of  the  total  bbni(650)  would  be 
associated  with  submicrometer  particles,  compared 
to  —45%  for  total  6,„(650).  Given  the  strong  depen¬ 
dency  of  bbp  on  composition  [e.g.,  Fig.  9(b)],  SAX  offers 
opportunities  to  advance  future  efforts  to  demon¬ 
strate  closure  in  forward  method  for  estimates  of  bm 
and  bhm  with  bulk  measurements,  particularly  for 
systems  enriched  with  minerogenic  particles. 

The  particulate  backscattering  ratio  ( bbp )  is  concen¬ 
tration  independent,  and  relevant  to  remote  sensing, 
light  field  modeling,  and  inference  of  bulk  refractive 
index  from  optical  measurements  [39,61].  Whitmire 
et  al.  [61]  presented  20%  as  a  conservative  estimate  of 
the  likely  maximal  error  of  bhp  based  on  a  propagation 
of  error  analysis  for  similar  field  instrumentation 
measurements  of  bulk  conditions  (absorption  and  at¬ 
tenuation  meter  over  nine  wavelengths,  ac-9,  Hy- 
droScat).  The  ratio  value  is  sensitive  to  composition 
(generally  higher  where  inorganic  particles  make 
greater  relative  contributions)  and  PSDs  [39,61]. 
While  the  data  set  presented  here  is  modest  in  size 
[Fig.  6(c)]  compared  to  recently  published  compila¬ 
tions  from  extensive  marine  surveys  of  the  ratio 
[39,61],  this  study  of  freshwater  systems  has  the  ad¬ 
vantage  of  more  definitive  specification  of  the  min¬ 
erogenic  component  of  bhp.  The  estimates  of  the  ratio 
for  the  minerogenic  particles  [ bbtfn ;  Fig.  6(c)]  represent 
upper  bound  values  of  bhp  corresponding  to  the  sce¬ 
nario  of  no  organic  (e.g.,  phytoplankton)  particles. 
The  modest  differences  in  the  bb  m  values  among  sites 
[Fig.  6(c)]  reflect  differences  in  PSDs  and  composition 
of  minerogenic  particles.  The  observed  6^(650)  val¬ 
ues  from  bulk  measurements  can  be  viewed  as  the 
outcome  of  the  mixture  of  6ftm(650)  values  with  the 
lower  values  associated  with  low  n  organic  particles. 
The  differences  in  the  bbp  values  were  generally  con¬ 
sistent  with  the  simple  view  of  the  mixing  of  different 
relative  amounts  of  minerogenic  and  organic  parti¬ 
cles,  as  indicated  by  the  dependency  of  bbp  on  the 
FSScorr/TSS  ratio  [Fig.  7(i)].  Loisel  et  al.  [39]  also 


reported  lower  bhp  values  for  particle  populations 
dominated  by  (organic)  material  of  low  refractive  in¬ 
dexes  and  high  values  for  sites  with  relatively  high 
concentrations  of  inorganic  particles.  The  population 
of  66p(650)  values  from  this  study  was  consistent  with 
those  reported  in  the  marine  compilations  [39,61] 
from  widely  different  scattering  conditions  (contribu¬ 
tions  of  inorganic  versus  organic  particles).  For  ex¬ 
ample,  our  mean  bhp  value  was  0.0132  compared  to 
0.0138  reported  by  Loisel  et  al.  [391,  and  a  geometric 
mean  of  0.013  reported  by  Whitmire  et  al.  [61].  More¬ 
over,  our  higher  values  were  within  the  upper  20%  of 
bbp  populations  in  the  marine  surveys  where  the  ele¬ 
vated  backscattering  ratios  were  attributed  to  enrich¬ 
ment  with  mineral  particles. 

E.  Advancing  the  Forward  Method  for  Light  Scattering 

Various  initiatives  would  improve  the  credibility  of 
estimates  of  the  minerogenic  component  of  light  scat¬ 
tering  based  on  SAX  data,  and  support  the  approach 
to  optical  closure.  Paired  measurements  of  (3P(B,  X) 
[62]  would  more  rigorously  define  bp(\)  and  particu¬ 
larly  fe6p(X),  and  the  limitations  of  Mie  theory  esti¬ 
mates.  Such  information  would  be  especially  valuable 
for  cases  where  minerogenic  particles  are  dominant. 
The  combination  of  well-defined  particle  morphome¬ 
try  from  SAX  and  detailed  P;,(0,  X)  information  could 
guide  the  development  of  empirical  adjustment 
approaches,  or  application  of  more  complex  formu¬ 
lations  of  light  scattering  [50],  for  cases  where  ap¬ 
plication  of  Mie  theory  to  estimate  scattering  is 
inappropriate  or  found  to  be  flawed.  Improved  quan¬ 
tification  of  PAV„,  and  thereby  60(X),  through  modifi¬ 
cation  of  SAX  techniques  would  improve  closure  for 
overall  bp(\).  For  example,  heavy  metal  staining  has 
been  used  to  improve  detection  and  resolution  of  or¬ 
ganic  particle  features  with  SAX  [63].  Peng  and  Effler 
[15]  have  described  approaches  to  reduce  the  modest 
uncertainty  associated  with  the  lying  flat  effect  on 
SAX  analyses. 

To  assess  the  applicability  of  the  findings  reported 
here,  minerogenic  PSDs  that  extend  into  the  submi¬ 
cron  range  should  be  documented  with  SAX,  and  Mie 
theory  estimates  should  be  made,  in  concert  with 
bulk  scattering  measurements,  for  a  wide  range  of 
freshwater  and  marine  systems.  Investigations  focus¬ 
ing  on  conspicuous  light-scattering  signatures  im¬ 
parted  by  minerogenic  particles  from  runoff  (clay 
minerals)  or  “whiting”  (CaC03  precipitation)  events 
represent  opportunities  to  advance  testing  of  the  ap¬ 
proaches  presented  here. 
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